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ABSTRACT 

The accreting millisecond pulsar SAX J1808.4— 3658 has shown a peculiar orbital evolution in the 
past with an orbital expansion much faster than expected from standard binary evolutionary scenarios. 
Previous limits on the pulsar spin frequency derivative during transient accretion outbursts were 
smaller than predicted by standard magnetic accretion torque theory, while the spin evolution between 
outbursts was consistent with magnetic dipole spin-down. In this paper we present the results of a 
coherent timing analysis of the 2011 outburst observed by the Rossi X-ray Timing Explorer and extend 
our previous long-term measurements of the orbital and spin evolution over a baseline of thirteen years. 
We find that the expansion of the 2 hr orbit is accelerating at a rate ~ 1.6 x 10~ 20 ss~ 2 and we 
interpret this as the effect of short-term angular momentum exchange between the mass donor and 
the orbit. The gravitational quadrupole coupling due to variations in the oblateness of the companion 
can be a viable mechanism for explaining the observations. No significant spin frequency derivatives 
are detected during the 2011 outburst (|z>| <4x 10 _13 Hzs _1 ) and the long term spin-down remains 
stable over thirteen years with v ~ — 10~ 15 Hzs -1 . 

Subject headings: binaries: general — stars: individual (SAX J1808.4— 3658) — stars: neutron - 
stars: rotation — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

The transient X-ray binary SAX J1808.4-3658 is the 
first acc reting millisecond X-ray pulsa r (AMXP) dis- 
covered (| Wijnands fc van der Klisl 11998) among the 14 
systems currently known. It is also the best sampled 
AMXP thanks to its relatively short recurrence time (1.6- 
3.3 yr) and the continuous coverage of the Rossi X-ray 
Timing Explorer (RXTE ) which has extensively mon- 
itored all outbursts since 1998. The presence of pulsa- 
tions reveal the spin of the accreting pulsar (« 401 Hz) 
and allow the study of torques that act upon the neu- 
tron star. This is particularly valuable because the 
spin evolution reveals details of the recycling mechanism 
that transforms a newly born slowly rotat ing neutron 
star into a millisecond pulsar via accreti on (|Alpar et al.l 
119821: iRadhakrishnan fc Srinivasa nlll982l). In our previ- 
ous wor ks (jHartman et al.1 120081: lHartman et alj l2009h 
we reported upper limits on the spin-up of the pulsar 
due to accretion torques (during outbursts) with a pul- 
sar spin frequency derivative of |z>| < 2.5 x 10~ 14 Hzs -1 . 
This value is smaller tha n predicted by accretion the- 
ory ()Ghosh fc Lamblll979l) if the spin frequency of SAX 
J1808. 4-3658 substantially differs from the equilibrium 
spin frequency at the accretion rate close to the peak of 
the outburst. The long term spin evolution of the pulsar 
reveals a constant spin-down of magnitude —5.5 ± 1.2 x 
10 _16 Hzs _1 possibly due to magnetic-dipole radiation 
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acting during quiescence 6 , for a surface magnetic field of 
the pulsar B ~ 1.5 x 10 8 G, in line with the expected field 
strength of millisecond radio pulsars. 

SAX J1808.4— 3658 is undergoing an unexpectedly 
fast orbital evolution with the orbi tal period in- 
creasing on a timescale of ~ 70Myr ( Hartma n et al.1 
I2008lldi Salvo et al.1 120081). The binary has an or- 
bital period of 2.01 hr (|Chakrabartv fc Morgan! 119981) 
and the donor star is a 0.05 — 0.1 Mp, brown dwarf 
(jBildsten fc Chakrabartvl 120011: iDelove et al.l I2008Q sug- 
gesting that the orbital evolution should be dominated 
by angular momentum loss vi a gravitational waves and 
possi bly by magnetic braking ()Tauris fc van den Heuvel 
2006). The timescale of the orbital evolution is, how- 
ever, too fast to be explained with such a scenario, 
and non-conservative processes with la rge mass loss 
from the system h a ve b een invoked (Idi Salvo et al.1 
20081 iBurderi et al.1 [20091) . lHartman et Sr HQOl j 
Hartman et al.1 ( 2009(1 suggested instead that inter- 
changes of angular momentum between the compan- 
ion and the orbit can dominate the short-term orbital 
evolution as seen in seve r al binary millise cond pulsars 
(|Arzoumanian et al.l 119941; iNice et al.ll2000D . 

On October 31 2011 Swift-B AT detected a new out- 
burst of SAX J18 08.4-3658 ([Markwardt et al.l [20Tlt 
iPapitto et al.l 120111) This is the 7th outb urst observed 
since its discovery (|in 't Zand et al.1 Il998|) and the 6th 
monitored with RXTE. We present a coherent pulsation 
analysis of the outburst and we complete the study of the 
spin and orbital evolution of SAX J1808.4— 3658 over a 
baseline of thirteen yr. 

2. X-RAY OBSERVATIONS AND COHERENT ANALYSIS 

6 In the abstract of Hartman et al. 2009 a wrong value of — 5.5± 
1.2 X 10 -18 Hzs - 1 is reported, due to a typographical error. 
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We use all RXTE Proportional Counter Array (PCA; 
I Jahoda et al.1 12006) public data for the 2011 outburst 
(Program-Id 96027). We construct the 2-16 keV X-ray 
light-curve with PCA Standard2 data averaging the flux 
for each observa tion and normalizing it in Crab units 
(see for example Ivan Straaten et ail 12003). One burst is 
detected at MJD 55873.9 and all data with a flux more 
than twice the pre burst level are removed from the light- 
curve. 

For the timing analysis we use all photons (excluding 
the burst interval) in the energy band « 2 — 16 keV (5-37 
absolute channels) in Event 122/is mode. The data are 
barycentered with the FTOOL fax bary by using the op- 
tical position (jHartman et al.ll2008t) and the JPL DE405 
solar system ephemeris. We fold ~ 500s long data seg- 
ments in pulse profiles of 32 bins, keeping only those with 
signal-to- noise > 3.3er, giving < 1 false pulse detection 
for the entire outburst. The S/N is defined as the ratio 
between the pulse amplitude and its la statistical error. 
The foldi ng procedure u s es the preliminary ephemeris re- 
ported in iPapitto et al.l (|2011h - A fundamental {v) and 
a first overtone (2v) are detected in the pulse profiles. 
The TOAs are measured separately for each harmonic 
to avoid that pulse shape variabili ty affects the fiducial 
point defining the pulse TOA (see iHartman et~aT1 [20081 
for details). 

To follow the evolution of the orbit and the pulsar spin 
we fit separately the two sets of TOAs (fun damental and 
first o vertone) with the software TEMP02 (jHobbs et al.1 
2006) . The initial model used is a Keplerian circular or- 
bit and a constant pulse frequency. We then repeat the 
folding procedure with the new timing solution until we 
reach convergence to the final orbital and pulse parame- 
ters. To calculate the errors on the rotational parameters 
we use Monte Carlo (MC) simulations, that account for 
the pre sence of long-timesca l e correlations in the puls e 
TOAs (jHartman et all 120081 : [Arzoumanian et al.lll994j ). 
To verify the presence of a spin frequency derivative 
we fit a pulse frequency and its time derivative to the 
TOAs and run 10 4 MC simulation s to estimate the sig- 
nific ance of the measurem ents (see IHartman et aLll2008l 
and lPatruno et all ((2009a) for further details). 

3. RESULTS OF THE 2011 OUTBURST 

3.1. X-ray Light-curve and Pulse Profiles 

The first RXTE pointed observation was taken on 
November 4, pa 5 days af ter the beginning of the outburst 
(|Markwardt et al"1l2011| ). The X-ray flux shows a peak 
at w 80 mCrab, remarkably higher than the peak lumi- 
nosity i n 2005 and 2008 and similar to the 1998 and 2002 
values (| Wiinandsl 12004) . Since the observations started 
5 days after the onset of the outburst, the true outburst 
peak at MJD 55868 (observed by Swift-BAT) has been 
missed by RXTE . The RXTE PCA light-curve of the 
2011 outburst is shown in Figure [TJ 

The flux is observed to slowly decay over the entire 
duration of the main outburst, reaching a minimum flux 
of 8 mCrab on MJD 55880. The flux then keeps de- 
creasing with two bumps at MJD 55881 and 55884, 
which might possibly be associated with the flaring- 
tail stage, when quasi periodic bumps are observed 
at low flux level j Wiinands et al.l 120011 : IHartman et al.l 
I2008t iPatruno et al.ll2009bD . However, differently from 
the 2000 to 2008 outbursts, this phase is poorly sam- 
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Fig. 1.— Top panel: X-ray light-curve (2-16 keV) of the 2011 
outburst (solid line) compared to the 2008 outburst (dotted curve). 
The data points are RXTE observations-long averages. Middle 
Panel: timing residuals for fundamental (blue dots) and first over- 
tone (red squares) for a v = const model. Each data point is an 
RXTE orbit-long average. The fundamental has a jump of ~ 0.15 
at MJD ss 55875. Bottom Panel: Fractional amplitude of the 
fundamental (blue dots) and first overtone (red squares). The 95% 
c.l. upper limits are indicated with pink circles (fundamental) and 
black open squares (first overtone). 

pled and difficult to characterize. During the tail 
a strong 1 Hz modulation was reported for several 
outbursts (2000, 200 2, 2005; Ivan der Klis et al.l 120001: 
IPatruno et al.l l2009b|) but it is not detected in any of 
the 2011 observations. 

The source became undetectable on MJD 55885.9 
(November 20) and the monitoring abruptly ended on 
November 27, due to solar constraints. These also pre- 
vented to establish whether at the time of the last obser- 
vation SAX J1808.4— 3658 was in quiescence or in one of 
the faint (~ 10 32 ergs _1 ) states observe d in the past with 
Swift and XMM-Ne wton observations (| Wi j nandsl 120031 : 
ICampana et al.ll2008l ). 

The pulse profiles are very sinusoidal until MJD 55874 
and then change becoming skewed with a more evident 
first overtone peaking on the right part of the profile. At 
MJD 55880, during the possible flaring stage, the pulse 
profiles become very sinusoidal again, a phenomenon 
never observed in the previous outbursts. This is re- 
flected in the rms amplitude of the two harmonics (see 

FigHD. 

3.2. Timing Noise and Error Estimation 

The timing solution of the 2011 outburst is displayed in 
Table [1] while a previous analysis of the five outbu r sts ob - 
served with RXTE is available in IHartman et al.1 (|2008h : 
IHartman et al.1 (|2009t ) . The authors found strong timing 
noise operating on the same timescales over which the 
pulse frequency and its time derivative were measured. 
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TABLE 1 

SAX J1808. 4-3658 Timing Solution for the 2011 Outburst 



Parameter 



Value 



Stat. Error Syst. Error 



v [Hz] 400.97520981 7 x 10" 8 10" 7 

HIHzs- 1 ] < 4 x 10- 13 (95% c.l.) 

P b [s] 7249.162 0.003 

A x sin(i) [lt-ms] 62.798 0.005 

T a sc [MJD] 55896.895635 0.00002 

e < IP' 4 (95% c.l.) 

The 2011 TO A residuals of a v = const model, show the 
typical behavior observed in the previous outbursts in 
both harmonics. 

We detect a strong phase jump in the fundamen- 
tal with magnitude of 0.15 cycles (0.3 ms) at MJD 
?» 55874 — 55876. This behaviour is similar to what was 
previously seen during the 2002 and 2005 outbursts, with 
phase jumps of 0.2 cycles observed when t he flux reached 
the transition from slow t o fast decay jBurderi et"al"1 
120061 lHartman et all 120081 iPatruno et all l2009bl) . The 
first overtone instead has no phase jump but displays 
a short timescale (few minutes to ~ 1 day) scattering 
slightly in excess of that expected from measurement er- 
rors alone. 

We use the phase information of the first overtone to 
pha se connect across the p hase jump of the fundamental 
(see lHartman et al.l 120081 for details of the procedure) . 
A net spin frequency derivative (i.e., measured over the 
entire outburst length) is not detected, with upper limits 
of \i>\ < 8.8 x lO^Hzs- 1 at the 95% confidence level. 
By removing the TOAs of the first observation in both 
harmonics (when timing noise is strong) and the 0.15 
cycle phase jump for the fundamental, we obtain simi- 
lar results with a more stringent constraint on the spin 
frequency derivative: \u\ < 4 x 10 _13 Hzs _1 (95% c.l). 

The errors on the orbital parameters measured with 
the fundamental are only marginally affected by timing 
noise since they are measured on timescalcs (2 hr) differ- 
ent than the timing noise one (~ days). This is verified by 
calculating a power spectrum of the TOA residuals and 
comparing the Poissonian level to the power at the or- 
bital frequency 1/P&. The excess power at 1/Pb is about 
1.5 and 3 times the Poissonian level for fundamental and 
first overtone respectively. We therefore rescale the sta- 
tistical errors on the orbital parameters by the same fac- 
tor. 

3.3. X-ray Flux - Pulse Phase Correlation 

In 2009 (|Patruno et al.l [2009c) proposed an alterna- 
tive method to partially account for the timing noise 
in the TOAs of AMXPs. X-ray flux variations were 
found to be linearly correlated or anti-correlated with 
the pulse phases. Instead of minimizi ng the rms of 
the TOA residuals. IPatruno et all (|2009c[) minimized the 
X 2 of a linear fit to the phase-flux correlation, finding 
slightly different spin frequencies than those measured 
with rms minimization methods. The reason of this 
difference is that instead of treating timing noise as a 
red noise process of unknown origin, the variations of 
the X-ray flux are assumed to instantaneously affect the 
pulse phases. We found that while the fundamental fre- 
quency always follows a correlation, the first overtone in 
some cases behaves differently and we exclude it from 
our 2011 analysis. By repeating the same procedure out- 



lined in IPatruno et all (|2009cT> for the 2011 outburst, we 
find a pulse frequency of v = 400.97520981(7) Hz, where 
la errors h ave been rescaled by a factor 2.5 such that 
X 2 /dof = 1 ([Bevington fc Robinsonll2003l : IPatruno et all 
2009d). The difference between this value and the pulse 
frequency v rms obtained with standard rms minimization 
of the TOA residuals is v — v rms = —0.15 ± 0.08/iHz. 

4. RESULTS ON THE LONG TERM EVOLUTION OF SAX 
J1808.4-3658 

4.1. Long Term Spin Frequency Evolution 

We first fit the change of the six constant pulse fre- 
quencies (from 1998 to 2 011) of SAX J1808.4- 3658 us- 
ing the values reported in IPatruno et all (|2009d ) and the 
2011 frequency obtained from the flux-phase correlation 
technique. We rescale the errors of each spin frequency 
to give a % 2 /dof = 1 and we fit a linear relation to the 
data. The fit gives a x 2 = 5.4 for 4 dof, and a spin-down 
of v = -1.65(20) xl0~ 15 Hz s" 1 (Fig. H). This is i n agree - 
ment with the value reported in IPatruno et all (|2009c). 
A spin frequency second derivative is not required by the 
fit, and we can place upper limits of \v\ < 10 _24 Hzs -2 
(95% c.l). 

We also fit the 2011 outburst pulse frequency with 
the rms minimization method (i.e., TEMP02 plus 
MC errors), together with the previous five spin fre- 
quency measurements obtained with the same technique 
(|Hartman et all 12009ft . The X 2 is high, 19.5 for 4 dof 
and we ascribe this almost exclusively to t h e 200 out- 
burst spin frequency (see lHartman et all (|2009h ). If 
we remove the 2000 data, the fit returns a x 2 — 3.57 
for 3 dof. The spin-down is constrained to be v = 
-7.4(4) x 10 _16 Hzs _1 . This esti mate is within 2a from 
the value previously reported bv lHartman et all (|2009f ) 
and within la when removing the 2000 outburst data. 

Since both fits of the long-term spin frequency evolu- 
tion are statitically acceptable, we cannot decide which 
of the two values reported is closer to the true v. The 
differences in the v and long term v found with the rms 
minimization and with the phase-flux correlation reflect a 
systematic uncertainty of ~ 0.1/i Hz that needs to be con- 
sidered until the exact mechanism behind timing noise is 
identified. 

4.2. Orbital Evolution 

To detect variations of the orbital period we us e 
the procedure expla ined in lHartman et all (|2008f k 
lHartman et all (|2009f) that requires estimates for the 
times of passage through the ascending node T asc at vari- 
ous outbursts. We c hoose as a refere n ce poi nt the T asCire f 
value in Table 1 of lHartman et all <|2009f ) . and we cal- 
culate the residuals AT asc = T asCii - (T asc , re / + N P b ), 
where T OSC) j is the i — th outburst and N is the closest 
integer to (T asCii — T„ „ r . re j) / Ph. The chose n Pi, is also 
reported in Table 1 of lHartman et all (12006" 

l200l 



In a previous work 

lHartman etall l2009f ) 

• m-i: 



we 
,-1 



(|Hartman et all 
detected a Pb = 
see also Idi Salvo et all 
When including the 2011 



(3.80 ± 0.06) x 10 
120081 and lFurderi et alj[2009h 

data, a parabolic fit gives a x 2 = 62.8 for 3 dof, thus 
suggesting that a constant increase of the orbital period 
is not the correct model for SAX J1808.4-3658 . We 

then add an orbital period second derivative (PfeJ and 
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1998-2011 Orbital Solution 




1998.3 2000.1 2002.8 2005.4 2008.7 2011.8 

Outburst Year 



Fig. 2. — Orbital (top panel) and spin frequency evolution (bot- 
tom panel) over 13 years of observations. The T a3c residuals cannot 
be fitted with a parabola (dotted line) and they need a cubic fit sug- 
gesting an acceleration of the orbital period = 1.65x 10 -20 s s — 2 . 
The spin-down of the pulsar is stable and proceeds at a rate of 
-10~ l5 Hzs" 1 . The reference frequency is u = 400.975210 Hz. 

fit a cubic polynomial to the data. The fit is statistically 
acceptable, with \ 2 — 5.4 for 2 dof and parameters 
P b = 3.5(2) x lCT^ss- 1 and P b = 1.65(35) x lCT 20 ss- 2 
(see Figure [2]) . 

A third time derivative or a sinusoidal model are not 
required by the data. We also tried to fix the P b and P b at 
the values observed up to 2008 and add a sinusoidal fit to 
the data to investigate the possibility of apsidal motion. 
The fit is statistically unacceptable and the periodicity 
found is of the order of 10 3 yr, which is already orders of 
magnitude larger than the expected general relativistic 
effect alone. 

5. DISCUSSION 
5.1. Pulsar Spin Evolution 

The long term spin-down continues with a constant 
rate comparable to what has been measured between 
1998 and 2008. The only plausible explanation for 
the spin-down in SAX J1808.4— 3658 requires magnetic- 
dipole radiation with a surface magnetic field of the neu- 
tron star at the poles B « (1.5 — 2.5) x 10 8 G (for a 
radius R — 10 km) for a magnetic-dipole moment fj, = 
(0.7 — 1.5) x 10 26 G cm 3 . This range includes the different 
spin-down measurement methods reported in § 13.21 and 
§ !3.3l and is close to the B field obtained with Fe line spec- 
tral fitting (jCackett et al.ll2009t IFapitto et aJj|2009ft an d 
accretion disk modeling (jlbragimov fc Pout anenl [20091. 

The remarkably constant long term spin-down places 
stringent constraints on any ongoing spin-up during an 
outburst. If spin frequenc y variations were large r than 
the upper limit reported bv lHartman etaD(|2Q08|) (\u\ < 
2.5 x 10~ 14 Hz s -1 ), they would produce a scatter in the 
observed spin frequencies of the order of 0.1 — 0.4/iHz (for 



v ~ 5 — 10 x 10~ 14 Hz s -1 ). This scatter is not observed 
with statistical errors of k, 0.05 — 0.1/uHz (see Fig [2]), 
and suggests that the net spin-up during an outburst is 
at best very small in magnitude for SAX J1808.4— 3658 . 

5.2. Orbital Period Evolution 

The orbital period P b indicates that the orbit of 
SAX J1808.4— 3658 is expanding and accelerating at a 
very fast rate. The timescale for the acceleration is: 

~ 2pb in m 
Tacc = -p- ~ wyr (1) 

Pb 

suggesting that the identified P b might not represent 
the secular evolution of the orbital period. If we as- 
sume that the measured acceleration is constant, then 
SAX J1808.4-3658 has changed sign of P b about 25 
years ago. There is, however, no reason for P b to be con- 
stant, since we are not sensitive to higher order deriva- 
tives and both P b and P b might be part of long timescale 
variatio ns similar to thos e observed in binary millisecond 
pulsars (jNice et alJl2000t ) . Until the 2008 outburst, when 
only a constant P b was detected, two interpretations were 
given: a secula r orbital evolution due to non-conservative 
mass transfer (jdi Salvo et al.ll2008l iBurderi et alJl2009t ) 
and a short-term evolution associated with exchange of 
angular momentum between the donor star and the orbit 
(|Hartman et alj|2008t [Hartman et al.ll2009D . 

If SAX J1808.4— 3658 has increased its mass loss due 
to an enhancement of the donor ablation, then the wind 
loss from the companion (M w ) needs to increase at a rate 
comparable with that of the orbital period, since P b oc 
M w . The mass loss is related to E ab i — 0.25 (R/A) 2 E, 
where E ab i is the ablation power, E is the pulsar rota- 
tional spin-down power and R and A the donor radius 
and the semi-major axis of the binary. To explain the 
acceleration P b , the energy loss of the pulsar E needs to 
have increased (in absolute value) in the last 13 years by 
a factor of w 5. Since E ex vu, the spin-down v needs 
to vary at a rate v ~ — 10 _23 Hzs -2 to reach the energy 
loss required. By using the long-term spin evolution pre- 
sented in the previous section we can put constraints on 
\i>\ < 10 _24 Hzs -2 . We conclude that the enhanced ab- 
lation scenario is not supported by the observations. 

A dynamically induced period derivative in the grav- 
itational potential well of a third body can also be ex- 
cluded. The effect of a potential well is identical on the 
orbital and spin frequencies and derivatives: 

/<*) = E (2) 

c 

where is the n-th time derivative of the orbital or 
spin frequency, a is the acceleration due to the third 
body, n is a unit vector along the line of sight and c 
the speed of light. To explain the observed P b we need 
a - 10" 15 - 10- 16 ms- 3 and ~ 10~ 21 Hzs-\ which 
is not observed. 

If the measured orbital evolution is a short-term event, 
then one explanation can be found with the donor spin- 
orbit coupling model. A coupling between the pul- 
sar rotational energy loss (in form of winds or fields) 
and the orbital angular momentum (|Damour fc Taylor! 
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1991) is ruled out by the small magnitude of the ef- 
fect produced by the tiny E of SAX J1808.4— 3658 . 
A mass quadrupole variation of the donor star is a 
more promising possibility. A change AQ in the 
mass quadrupole leads to a change in orbital pe- 
riod (iRichman et al.1 Il994t LApplegate fc Shahaml 11994 
lAppleeate fc Pattersonlll987f) : 



AP, 



n 2 R 3 m s An 

GM ~M~~tt~ 



(3) 



where M and R are the donor mass and radius, M s is 
a thin shell of mass generating the quadrupole, and fl 
the angular velocity of the star. If we assume that the 
angular velocity of the donor is almost synchronous with 
the orbital angular velocity, then the variation AP,, ~ 
0.004 s observed in the last 13 years gives: 



Afl 



t M 
M7 



(4) 



The observed orbital period variations in the eclipsing 
millisecond pulsar PSR J2051-0827 and PSR B1957+20 
are likely to be caused by changes in the quadrupole 
moment of the companion (jArzoumanian et al.1 



1994t iDoroshenko et all I200H lLazaridis et"aH 120111) . 
ApplegateT i 19921 ) proposed a magnetic activity cycle 
that leads to a deformation of the star at the origin of 
this behavior. The donor star of SAX J1808.4-3658 is 
also in Roche lobe contact, whereas binary millisecond 
pulsars are detached systems. If the orbital period 
of SAX J1808.4— 3658 has decreased in the past for 
some time, then the Roche lobe has moved across the 
outer envelope of the brown dwarf enhancing the mass 
transfer rate. A detailed discussion of this effect is be- 
yond the scope of this letter, but we can speculate that 
SAX J1808.4— 3658 has gone through periodic episodes 
(each lasting r acc ~ 10 yr) of enhanced accretion in 
the past. This effect is opposite during the accelerated 
orbital expansion, with the mass transfer being less than 
in the non-accelerated case. The quadrupolar moment 
change has also the effect of heating the star, providing 
an ex planation for the large entropy content of the 
donor (pelove et al.ll200l) . 
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